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ABSTRACT 

Voltage-driven spin transfer torques in magnetic tunnel junctions provide an outstanding tool to design advanced 
spin-based devices for memory and reprogrammable logic applications. The non-linear voltage dependence of the 
torque has a direct impact on current- driven magnetization dynamics and on devices performances. After a brief 
overview of the progress made to date in the theoretical description of the spin torque in tunnel junctions, I present 
different ways to alter and control the bias dependence of both components of the spin torque. Engineering the 
junction (barrier and electrodes) structural asymmetries or controlling the spin accumulation profile in the free 
layer offer promising tools to design efficient spin devices. 
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1. INTRODUCTION 

While most of the commercial microelectronic devices are based on the charge of the carriers (electrons and holes), 
Spintronics relies on the spin angular momentum of the carrier to generate low energy consumption functional 
devices. The most impressive demonstration of the technological relevance of Spintronics is the implementation 
of Giant Magnetoresistive read heads in magnetic data storage.^ As illustrated by the current state of the art 
of data storage, the first virtue of magnetic devices is their non- volatility: in the absence of external input 
and thermal activation, the magnetic state remains frozen in its equilibrium direction. Therefore, non-volatile 
Magnetic Random Access Memories (MRAM) composed of magnetic tunnel junctions are expected to present 
outstanding performances down to 30nm and below. 

There are currently various types of MRAM. The original MRAM concept, commercially available, is based 
on field-induced switching but has limited performances when scaling down the device areaP Another candidate 
is the thermally-assisted MRAM, or TAS-RAM, where field-induced switching is assisted by thermal activation. 3 
Beyond these different alternatives, spin transfer-MRAMs (STT-MRAMs) present a particular interest for both 
applied and fundamental condensed matter physicists^ through the electrical manipulation of the magnetic 
configuration of the device. This has been made possible by the prediction of the so-called spin transfer torque 
(STT) by Slonczewski and BergerP A spin-polarized current impinging on a ferromagnet can transfer its spin 
degree of freedom to the local magnetization exerting a torque on it. This effect has been observed in metallic 
spin- valves (two ferromagnetic layers separated by a spacer), tunneling junctions and even magnetic domain 
walls and is shown to induce current- control magnetic excitations, GHz self-sustained magnetic precessions and 
magnetization switchingPE' This area has known tremendous developments in the past ten years with the 
realization of efficient devices such as spin torque MRAMs, radio- frequency oscillators and reprogrammable 
logics. 

Beyond technological challenges, spin transfer torque in magnetic tunnel junctions (MTJs) presents a number 
of puzzles. Although STT in MTJs is considered to have reached maturity, the bias dependence extracted from 
routine measurements is still not well understood nor controlled. In this article, I propose an overview of the 
means to control the bias dependence of STT in MTJs. I present first the nature of the spin torque expected in 
an "ideal" ballistic tunnel junction. Then, I briefly review the influence of asymmetries and interfacial scattering 
and finally I discuss the influence of spin diffusion in the free layer. 
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2. STATE OF THE ART: EXPERIMENTS 



Spin transfer torque in tunnel junctions has been first demonstrated in MgO-based junctions by Sun and in AlOx- 
based junctions by Huai et al. and Fuchs et al. 8 However, Fe/MgO/Fe-type junctions very rapidly supplanted the 
AlOx junctions due to their outstanding transport properties. Since then, a number of exciting realizations have 
been achieved including switching, excitations, precessions, spin diode etc. We refer the reader to Ref. [6] for 
more details about the experimental achievements. Among the most important results, it has been demonstrated 
experimentally and theoreticalljDMH that the spin torque has the following form 



Here, m and p are the magnetization direction of the free and pinned layers, respectively and V is the bias 
voltage applied across the junction. The first term Ty in Eq. ([!]), referred to as the in-plane torque, competes 
with the magnetic damping allowing for self-sustained magnetic precessions and switching, whereas the second 
term Tj_, referred to as the out-of-plane torque, acts like an effective field applied along p. As we discuss in the 
next section, this dependence is theoretically predicted in ballistic junction when the transverse spin density is 
absorbed at the interface. 

However, recent experiments have reported important discrepancies between the actual bias dependence 
of the spin torque and the one proposed in Eqs. Q-Q. 17 20 In particular, Oh et a/P^l showed that in an 
asymmetrically designed MTJ, the bias dependence of the out-of-plane torque acquires a linear contribution. 
This linear contribution has also been observed by Petit et al. and Heinonen et alP^ as small biases. Deac et 
alP^ and Li et alP^ have also uncovered complex bias dependences at large bias voltages. The aim of the present 
article is to present a coherent description of the spin transfer torque in tunnel junctions in order to account for 
these different observations. 



At this early stage of the theory, we directly connect the spin torque acting on the volume of the free magnetic 
layer to the interfacial spin current transverse to the local magnetization, = Js^ {interface- By definition, 

the spin transfer torque is due to the amount of spin angular momentum transferred to the magnetic layer, 
minus the amount of spin lost through spin-flip mechanisms. In the present case of ballistic transport (no spin 
relaxation), incoming spins gradually align on the local magnetization by transferring its transverse momentum 
to it, as schematically depicted in Figjlja). This alignment takes place on a very small distance near the interface, 
named the spin dephasing length A^. Assuming that the spin torque is related to the interfacial spin current 
means that the thickness of the free layer is much larger than the spin dephasing length [d » A^). As we will 
discuss in the following, this definition is not appropriate in the general case and quantum resonances as well 
as spin diffusion effects can significantly affect the torque in the case of ultra thin free layers. Nonetheless, we 
assume for now that the spin transfer is interfacial and given by the interfacial tunneling spin current. 

The available theoretical descriptions of the bias dependence of spin torque in tunnel junctions have been ob- 
tained through transfer matrix formalism pi free electron JHHHHH1 tight bindin^ESEH and ab-initio calculations.^ 
Although these methods capture some or most of the band structure details of the system, they do not provide 
conceptual pictures about the symmetries of the bias dependence of the spin torque. In the present paragraph, 
we wish to provide such a qualitative description, disregarding any band structure details. A convenient method 
is the Bardeen Transfer Matrix (BTM) approach} 2 ^ that was widely used by Slonczewski in his descriptions of 
spin torque in tunnel junctions. 11 The tunneling transport is expressed as the product of the interfacial densities 
of states through a transfer matrix. The simplest version of this model proposed by Julliere 2 ^ had a reasonable 
success in explaining the tunneling magnetoresistance. Formally, BTM approach can actually be recovered as 
an asymptotic limit of the Keldysh tight-binding model in the case of thick and large barriers!^ As long as the 
interface states are decoupled (tunneling is treated perturbatively), BTM is applicable. Therefore, the present 
discussion does not reproduce quantitatively the physics of ultra thin barrier MTJs, but we believe it does provide 



T T\\m x (p x m) + T±m x p, 
T,| = ai V + a 2 V 2 , Tj_ = b + b 2 V 2 . 



(1) 
(2) 



3. BALLISTIC INTERFACIAL SPIN TORQUE 



3.1 Ballistic Tunneling 




Figure 1. (Color online) (a) Schematics of the spin transfer process in a magnetic tunnel junction. The black arrows refer 
to the magnetizations direction, and the red arrows represent the itinerant electron spins. Note that the itinerant spins 
get aligned on the local magnetization of the free layer m over a distance A^, called spin dephasing length; (b) Schematics 
of the alignment of the impinging spin in the rigid spin picture (red arrow) and in the realistic transport description (blue 
arrow). This supplementary angle gives rise to the out-of-plane torque. 



the relevant symmetries and tendencies. In the spin-dependent BTM formalism, the charge and spin current 
densities are given in the spinor form in the 2x2 spin space 

J = 2tt| J de\p L f L ^ R p R fl^ R f L (l - f R ) - PR f R ^ L p L f R _+ L f R (l - f L )], (3) 

where Jl,r and Pl,r are the Fermi distribution function and electronic density of states at the left (right) 
interfaces, and T^^r (Tr^l) is the spin-dependent transfer matrix accounting for both elastic and inelastic 
tunneling. In the spinor formalism, the charge current and spin current are expressed j e = Tr[J] and J s = 
Tr[<xJ], where & is the vector of Pauli spin matrices. The hat "denotes 2x2 matrices in spin space. This form 
accounts for contributions of rightward and leftward electrons originating from the left (L) or right (R) reservoirs. 
Using the above definition assumes that tunneling is only due to interfacial densities of states and transmission 
through the barrier. 

The important point here is the definition of the transmission matrices Tl^r and Tr^l. In the case of spin 
torque, we have to extend the BTM formalism to non-collinear magnetization directions. That is to say, the 
interfacial densities of states must be rewritten in the absolute quantization axis that we chose aligned on the 
magnetization of the right layer. Slonczewski^ and Levy and FertP^ used this approach to derive explicit low 
bias expressions of the in-plane torque. They assumed that the ballistic tunneling through the barrier is spin 
independent. In other words, the spin state does not rotate during the tunneling. The spin impinging on the free 
layer is then simply aligned on the polarizer orientation p, as depicted in Fig. [ljb) . This picture, that we refer 
to as the rigid spin assumption, allows to get an expression for the in-plane torque, but does not provide any 
indications about the nature of the out-of-plane torque. However, numerical calculation^EMSj and experimental 
evidencdlE3f!£] have demonstrated the existence of a large out-of-plane component of the spin torque. This 
indicates that a fundamental ingredient is missing in the original BTM approach and corrections to the rigid 
spin assumption must be considered. Actually, when a spin-polarized electron impinges on a ferromagnet, its 
transmission through the interface is accompanied by a spin rotation. This can be easily understood using a toy 
model. Let's consider a rightward electron in the rotating frame of the right electrode. Its original wave function 
|\I/)o is transformed into \^) t under transmission through the interface 

l*)o = cos^| t) -sin^| I) |*>, = f|*) = (Ttcos^| t) -^sin^| |». (4) 

Since the interface between the barrier and the free magnetic layer is spin-dependent, it filters the electron spin 
and its majority and minority projections are not transmitted the same way: the interaction with the interface 
induces a phase shift between the majority and minority projections of the spin. In other words, since the 
impinging spin feels a magnetic interaction at the interface, it rotates upon transmission/reflection. This effect 
has been identified by Stiles and Zangwill 26 in metallic spin- valves and by Manchon et al. 12 in magnetic tunnel 
junctions. These references demonstrate that the supplementary angle gained over transmission can be quite 



large and depends on the direction of incidence. In metallic systems, averaging over the Fermi surface quenches 
this contribution,^ whereas in magnetic tunnel junctions, the wave vector filtering yields an effective angle that 
is quite large and has a component out of the (m,p) plane^ [see Fig. [ljb)]. Consequently, the effective spin 
density in the free layer is rotated out of this plane and produces a so-called out-of-plane torque. 

This implies that the tunneling matrix is no more spin independent and that the incident spin "sees" a 
free layer magnetization with a virtual orientation ip being the phase induced by the spin-dependent 

transmission. The expression for the tunnel current spinor in the quantization axis of the right (free) layer is 
then 

J = 2tt| y de[TZp L f L ^ R TZ VL p R nl L fl^f L (l - f R ) - p R t R ^ L K VR p L K%Ji_+ L f B {l - f L )), (5) 
where we defined the interfacial density of state of the ith electrode and the rotations matrices 1Z and 



If pi + Api \ / cos 



up < 



Pi ~ 2\ Pi-Apt J ' K ~ \ sinf cosf J ' ^ " V e" 4 * sin § cos f J' W 

As mentioned above, note that the matrices R lfi account for the rotation of the electron spin when tunneling 
through the barrier whereas 7Z ensures that the spinor current is defined in the quantization frame of the right 
layer. A straightforward calculation leads to 

T N (e) = 2^Jde\T\ 2 p R Ap L (cos^ L f L -cos^ R f R ), (7) 
Tl(c) = 2^J de\T\\p R Ap L sm^ L f L + p L Ap R sm^ R f R ), (8) 

At this stage we have made no assumption about the symmetry of the junction. Note that in general the 
correction from the interfacial precession is small (p « 1. In this case, it clearly appears that the in-plane 
torque is antisymmetric (in the limit of small cp) when exchanging L and R indices, whereas the out-of-plane 
torque remains symmetric. We considered ballistic tunneling in a thick barrier, but it is in principle possible 
to have asymmetric electrodes or asymmetric barrier profile. Let's now assume a symmetric tunnel junction 
submitted to a bias voltage. In this case, we can rewrite the explicit energy dependence of the different terms 
involved in the spin transport 

Pi (e) = p(e ± eV/2), A Pl (e) = Ap(e ± eV/2), ^(e) = ^(e ± eV/2). (9) 

where the sign ± is associated with the left (right) electrode respectively, and the Fermi-Dirac distribution is 
given by /l,r(^) = (e^ e ~ €F±eV ^ 2 ^ kBT + l) -1 . Then, the in-plane and out-of-plane torques become 

T,| = J de[ttt(e + eV/2)-tU(e-eV/2)], (10) 

de[t ± {e + eV/2) + ^(e - eV/2)}. (11) 



Therefore, independently on the detail of the band structure and in the limit of the model, if the junction is 
symmetric, it implies that 

T|, = aiV, (12) 

T ± = b + b 2 V 2 , (13) 

in qualitative agreement with the numerical estimates E2EMH1 Note that the actual energy dependence of the 
interfacial spin-dependent densities of states influences the bias dependence of the spin torque at large voltages. 
First, as experimentally demonstrated by Valenzuela,^ the polarization of the impinging electrons depends on 
the bias polarity which induces quadratic corrections to the simple linear voltage dependence of the in-plane 
torque presented in Eq. [12] Second, at large bias dependence, Fowler-Nordheim processes start dominating the 
transport which induces oscillatory voltage dependences of the tunneling magnetoresistance and spin torque, as 
calculated by Tang et alEMSH 



3.2 Materials Consideration 



When a spin-polarized electron enters into the free layer, its spin precesses around the local magnetization with 
a spatial period of 2it/(k^ — fc^)HHH25l Since the angle between the electron spin and the local magnetization 
depends on the incident angle of the electron, the precession period also depends on the angle of incidence. After 
averaging over the Fermi surface, as mentioned above, the effective magnetization displays a oscillatory damped 
spatial profile as depicted schematically in Fig. [2] 

The most common system investigated to date is the MgO-based magnetic tunnel junctions with CoFeB 
electrodes. For thick enough barriers, the electron wave function is filtered so that the transport is dominated 
by Ai symmetries. 29 Since Ai bands are half metallic in CoFeB, the junctions display huge tunneling magne- 
toresistance ratios In this limit, the majority projection of the electron spin propagates regularly in the free 
layer whereas the minority spin wave function is evanescent. As a consequence, the spin density profile is heavily 
damped, as shown in Fig. |2j and the spin torque is very strongly localized at the interface. Thus, the model of 
interfacial spin torque proposed above is valid. Note however that in thin MgO barriers, resonant states arise 31 
that alter the tunnel magnetoresistance by allowing for minority electrons from other band symmetries to tunnel 
through the barrier. In this case, the minority waves are no more evanescent and are allowed to propagate in 
the free layer, yielding a torque that can extend through the volume of this layer. 




Figure 2. (Color online) Schematics of the spatial profile of the spin density transverse to the local magnetization, extending 
in the free layer. We represent typical profiles in the cases of a weak ferromagnet (long precession length - red), strong 
ferromagnet (short precession length - black) and half-metallic free layer (exponential decay at the interface - blue). 

3.3 Structural Asymmetries 

The previous simple model provides a qualitative description of the bias dependence expected in the case of 
asymmetric junctions. When the work function of the left and right electrodes are different or when the materials 
of the electrodes are different from each other, they induce asymmetries in the transport: the spin torque is 
more efficient at one bias polarity than at the other. In other words, when inserting structural asymmetries, 
tl^r 7^ tr^l- The first implication is that the bias dependence of the out-of-plane torque deviates from the 
"conventional" quadratic dependence. This is schematically represented in Fig. [3j left panel. In this figure, 
we represented the bias dependence obtained using a free electron model in the case of symmetric (solid lines) 
and asymmetric tunnel barrier (dashed lines). Note the presence of a quadratic component V 2 in the in-plane 
torque, that arises from the energy dependence of the interfacial densities of states. Details can be found 
elsewherelHEHlH 

This feature has been recently exploited by Oh et alP^l to reduce the back-hopping process observed at 
large biases.^ As discussed above, in a symmetric junction, whereas the sign of the in-plane torque depends 
on the polarity of current injection (ex V), the out-of-plane torque is always in the same direction (ex V 2 ). 
Consequently, at positive polarity both torques favor the antiparallel configuration, whereas at negative polarity, 
the in-plane torque favors the parallel configuration while the out-of-plane torque favors the antiparallel one. 
This competition leads to back-hopping of the magnetization state which is detrimental for applications such as 
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Figure 3. (Color online) (a) Bias dependence of the in-plane (black) and out-of-plane (red) torques in the case of symmetric 
(solid) and asymmetric barrier (dashed). See Refs. *** for details; (b) Spin torque magnitude as a function of the free 
layer thickness in the case of coherent ballistic transport (quantum interferences are visible - black), spin diffusion in the 
short spin dephasing limit (the thickness dependence in oc 1/d - red) and long spin dephasing (deviation from 1/d - blue). 



MRAMsllIll22l Introducing artificial asymmetries in the system, by using either different electrode materials or 
by engineering the barrier potential profile, Oh et al. 17 showed that it is possible to quench this back- hopping 
process by adding a linear part to the bias dependence of the out-of-plane torque. 

3.4 Interfacial Scattering 

We now briefly review previous works on the influence of interfacial spin scattering in magnetic tunnel j unc- 
tions E3H1 It is well known that electron interactions with thermal excitations can affect magnetic tunnel junc- 
tions properties such as the magnetoresistance!^ Hot electrons-induced magnon emission and absorption affect 
the polarization of the tunneling electrons and electron-phonon interactions result in thermally-assisted tunnel- 
ing. In the presence of electron-magnon and electron-phonon interactions, the transfer matrix becomes both spin 
and energy dependent 



f e L Z? R = f L ^ R [i+^(&.S«+<r.St r )\, (14) 




THr = r L ^[l+\/^_(6 q + 6j)|7, (15) 

where Q m (Q^ 1 ) is the phenomenological electron-magnon (electron-phonon) efficiency, N is the number of atoms 

per cell, a is the vector of Pauli spin matrices and S tr are the transverse part of the magnetizations of the 
left and right electrodes. Details about the derivation of Eqs. ([l4|)-([T5|) can be found in Ref. [25]. Assuming 
that the electron spin-dependent densities of state do not vary much over the range eV, and considering acoustic 
phonons (u oc q, Q q oc q) with a density of states of the form p p h(w) oc uj u , we obtain, at T=0 K and low bias 
voltage 

T N =G P L sin 9(l + ( ph | V\» +2 )V, (16) 
T± - T ±0 = G P R <pLsmO( ph \V\» +3 , (17) 

( p h being a coefficient that depends on the electron-phonon coupling, Fermi energy, Debye temperature etc... 
Note that the symmetry of the out-of-plane torque against the bias is conserved, whereas the in-plane torque 
acquires an antisymmetric component. 

In the case of electron-magnon interaction, the transfer matrix [Eq. ( [l4| )] possesses non-diagonal elements 
that are responsible for spin-flip scattering. We then expect a much more complex influence on the torque. 
Assuming a magnon density of states of the form p m (oj) = uo v ', symmetric electrodes {P l = P, <pi = cpn) and 



T=0 K, we find: 



T|| - T|| oc sin<9[P(l + P) - (1 - P)(l + Pcostf)] V^ 2 , (18) 
T ± - T ±0 oc P0sin<9(l - cos6>)|Fr +2 . (19) 

The detail of these expressions can be found in Ref. [25]. Interestingly the out-of-plane torque and the conduc- 
tance (not shown) both acquire a component that is symmetric against the bias because neither electron-magnon 
nor electron-phonon scattering break the junction symmetry. Furthermore, since the electron-magnon interac- 
tion mixes the majority and minority channels, the angular dependence is also affected, contrary to the case of 
electron-phonon coupling. 

3.5 Finite thicknesses in ballistic regime 

In the case where the electrodes' thickness is much larger than the spin dephasing length, the incident spin 
current is totally absorbed in the free layer. Therefore, the torque is determined by the interfacial spin current 
only. For example, in the case of half-metallic transport of Ai electrons in FeCo/MgO/FeCo junctions, since 
Ai minority electrons cannot propagate into the free layer, the length over which spin precesses (or in other 
words, the length over which spin transfer takes place) is reduced to one to two monolayers.^ In this case, the 
total torque exerted on the layer does not depend on the layer thickness since all the injected transverse spin 
current has been absorbed at the interface. Therefore, the average torque acting on the free layer is inversely 
proportional to the layer thickness (ex 1/d). 

In the case of normal metallic behavior, both majority and minority spin propagate in the free layer, giving 
rise to a non- vanishing coherence length (see Fig. 2). In the case of free layers of thicknesses comparable to the 
spin precession length, the transverse spin current is not fully absorbed in the free layer and is reflected by the 
second interface. This induces quantum interferences between the injected spin current and the reflected (counter 
propagating) one. The precession pattern of the spin density is therefore a combination between ex.p[i(k^ — k^)x] 
and exp[i(k^ + k^)x]. This interference pattern is therefore affected by the thickness of the free layer and 
produces a complex thickness dependence, as schematically depicted in Fig. |3jb). This dependence has been 
studied numerically in several references Ullll^l Note that quantum coherent is seminal to obtain this pattern and 
interfacial roughness or impurity scattering can easily destroy these interferences. 



4. DIFFUSION OF TUNNELING SPIN TORQUE 

In the present section, we briefly present recent results obtained in the opposite limit of diffusive transport in 
the free layerP^ As seen in the previous section, in the coherent ballistic regime, quantum interferences occur 
when the free layer thickness is finite. However, in the limit of diffusive transport, quantum interferences are 
quenched and spin relaxation and diffusion dominate the transport. The tunneling process through the insulating 
barrier imposes a ballistic injection of carriers at the interface between the tunnel barrier and the free layer which 
constitutes a boundary condition to the coupled diffusive spin transport equations for the transverse component 
of the spin accumulation vector s. Along these guidelines, the spin dynamics of the transverse spin accumulation 
in the free layer is governed by the following coupled equations in steady state 

V • J = — — s x m — —m x (s x m) — s, (20) 

70 v T s f 

J = -£>V(g)s, (21) 

where s is the spin accumulation, m is the direction of the local magnetization, and J is the spin current 
tensor. The diffusion is characterized by the diffusion constant V and the spin dynamics is controlled by the 
spin precession time tj, spin decoherence time and spin relaxation time r s f. Whereas the spin relaxation 
affects the three spin components, the spin precession and spin decoherence terms only affect the two transverse 
components of the spin accumulation vector. The spin torque is defined as the spatial change of spin current, 
compensated by the spin relaxation term 



i / dn(-V-J--s). (22) 



Here, Q is the volume of the magnetic layer. Note the seminal difference with the definition exploited in the 
previous section: the spatial variation of the spin current is now partially balanced by the loss of the spin angular 
momentum through spin flip scattering and the electron is now propagating diffusively in the free layer yielding 
a specific spatial distribution of the spin current that needs to be evaluated. 

We consider a finite free layer embedded between the tunnel barrier and a normal metallic capping layer. 
Assuming that a spin current Jo = J\\ + ij±_ is imposed by the ballistic tunnel at the interface with the barrier, 
and considering a standard continuity condition at the interface with the capping layer™ the total spin torque 
exerted on the ferromagnet is 



Jo L 2 cosh f 



cosh ^ 



r] sinh jr — 1 



- 77 sinh j£ 



(23) 



where ^2 = — \V + tt, A = — tV + A + t^-, d is the free layer thickness and ri is a parameter that only 

^0 A J A ^ A J A A sf 

depends on the bulk characteristics of the free and capping layers. We clearly see that in general the in-plane 
and out-of-plane torques are a mixture between in-plane and out-of-plane interfacial spin currents 



T\\=aJ\\+0J ± , 
T± = -PJ\\+aJ±, 



(24) 
(25) 



These expressions reduce to the ballistic case studied in the previous section in the case of infinitely free layer 
thickness (d — >• 00) and infinite spin diffusion length (X s f — >• 00). However, in the most general case, it appears 
that the spin torque can not be simply identified to the interfacial spin current. This property arises from the 
fact that (i) spin relaxation is always present, alters the propagation of the spin density in the metallic layers 
and redistributes the spin degree of freedom on the two torque components; (ii) when the free layer thickness 
is comparable to the spin dephasing length, the transverse spin density diffuses towards the interface with the 
capping layer which in turns influences the spin dynamics in the free layer. This last effect virtually enhances 
the impact of the spin diffusion length on the spin current mixing. These two situations are illustrated in Fig. [4j 
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Figure 4. (Color online) Spatial distribution of the out-of-plane (black) and in-plane (red) components of the spin accu- 
mulation in the free layer/capping layer bilayer structure for short spin dephasing length (A</> = 0.5nm - solid lines), long 
spin dephasing length (A</> = 1.5nm - dashed lines) for short spin relaxation length in the capping layer (A^ = lnm). The 
dotted lines display the spin accumulation profile for long spin relaxation length in the capping layer (A^ = 1.5nm and 
X^f = lOnm). The interface with the tunnel barrier is located at x — and the interface between the free layer and the 
capping layer is at x = 2nm. These curves were calculating following the theory developed in Ref. 34. 



The diffusion of spin current in the free layer has two major implications: the redistribution of the spin angular 
momentum is thickness-dependent and as soon as the thickness of the free layer exceeds the spin dephasing length, 
the interfacial spin current limit is recovered [see Fig. [3^b), red line]. Note that the deviation we obtain here 



is not related to quantum interferences (which are destroyed by impurity scattering and roughness), but rather 
to the incomplete absorption of the spin current: when the thickness of the free layer is smaller than the spin 
dephasing length A^, the transverse spin current responsible for the spin torque is not fully absorbed in the free 
layer and diffuses towards the capping layer. This induces a deviation from the usual 1/d-thickness dependence 
[see Fig. |3|b), blue line]. Increasing the thickness of the free layer improves the absorption of the spin current 
and for thicknesses much larger than the spin dephasing length, the thickness dependence of the torque recovers 
the 1/d limit. Note that in the case of half-metallic behavior, as in Fe/MgO/Fe tunnel junctions, the minority 
band of Ai symmetry does not propagate in the ferromagnet which induces a quenching of the spin dephasing 
length and reducing the spin torque to a 1/d behavior. 

Second, this new dynamics mixes the two transverse components of the spin current and is therefore respon- 
sible for the deviation from the ballistic bias dependence shown in Eq. ([2]). If one assumes a bias dependence 
of the spin current on the form J\\ — a±V and J±_ = 62 V^ 2 , as expected and observed for systems such as 
Fe/MgO/FeJSEH then both in-plane and out-of plane spin torque components will be a mixture of linear and 
quadratic bias dependences 

T,| = aaiV + /3a 2 V 2 , (26) 
T ± = -p ai V + aa 2 V 2 , (27) 

This spin mixing may explain the recently obtained experimental results that show a strong linear voltage 
dependence of the out-of-plane torque, in contrast with the predictions 

5. CONCLUSION 

In conclusion, the role of spin diffusion in the metallic layers of MTJs has been addressed theoretically. Assuming 
an interfacial bias-driven spin current at the interface between the insulator and the ferromagnet, the spin 
diffusion equation is solved and describes a complex spin dynamics in the metallic layers. It is found that 
this dynamics mixes the components of the spin current tranverse to the local magnetization which results in a 
superposition between linear and quadratic bias dependence of the out-of-plane torque. The thickness dependence 
of the spin transfer torque is also altered for small thicknesses. 
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